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The  multidimensional  solitons  in  a plasma:  structure,  stability  and  dynamics 

V.Yu.  Belashov18,  S.V.  Vladimirov2  b 
'Kazan  State  Power  Engineering  University,  51  Krasnosel'skaya,  Kazan  420066  Russia 
2School  of  Physics,  University  of  Sydney,  NSW  2006,  Australia 

The  formation,  structure,  stability  and  dynamics  of  multidimensional  nonlinear  waves  and  solitons 
in  a plasma  with  v = -(zx  V\jr)  / B and  P > 1 are  studied.  To  study  the  stability  of  multidimensional 
solitons,  the  variation  problem  of  the  Hamiltonian  bounding  with  respect  to  deformations  con- 
serving momentum  is  used.  To  study  evolution  of  solitons  and  their  collision  dynamics  the  equa- 
tions are  integrated  numerically.  It  was  obtained  that  in  both  cases  the  formation  of  multidimen- 
sional solitons  can  be  observed.  It  is  found  that  the  soliton  elastic  collisions  can  lead  to  formation 
of  complex  structures  including  the  multisoliton  bound  states. 


1.  Basic  equations 

In  this  paper,  we  study  formation,  structure,  stability 
and  dynamics  of  multidimensional  solitons  formed  on 
the  low-frequency  branch  of  oscillations  in  a plasma  for 
Ps4jw77B2  «land  P>1.  These  oscillations  are  de- 
scribed by  equation 

f X 

d,u  + A(t,u)u  = /,  / = K ALudx, 

~ ~ J-°°  U) 

A±  = dr  +dz. 

For  Aj(r,K)  = audx  -9x(v-P3x -yS3)  Eq.  (1)  falls 
into  GKP  (Generalized  Kadomtsev-Petviashvili)  class 
of  equations,  and  in  the  case  when  p = 4nnT  I B2  « 1 
for  co<(0#  = eB  / Me,  kXD  «1,  describe  propagation 
of  the  fast  magnetosonic  (FMS)  wave  in  a magnetized 
plasma  with  k2»k2,  vx«cA  near  the  cone  of 
0 = arctan  (M I m)U2  [1],  In  this  case,  the  function  u is 
the  dimensionless  amplitude  of  the  magnetic  field  of  the 
wave  h = B„/B,  the  factors  at  the  terms  describing 
nonlinearity,  dissipation  and  dispersion  effects,  respec- 
tively, are  defined  by  plasma  parameters  and  the  angle 

0 = (B,k) . In  opposite  case,  A2(t,u)  = 3s  I p\  udx- 

-d2(iX  + y/) , Eq.  (1)  converts  into  3D  derivative  non- 
linear Schrodinger  (3-DNLS)  equation  class  and  in  the 
case  when  p > 1 describes  dynamics  of  the  finite-am- 
plitude Alfven  waves  propagating  nearly  parallel  to  B 
for  u = h = (By+iBz)/2B\l-$\,  h = Bx/B0  where 

p = (1  + ie) , and  e is  the  “eccentricity”  of  the  polarizati- 
on ellipse  of  the  Alfven  wave  [2].  The  upper  and  lower 
signs  of  X = + 1 correspond  to  the  right  and  left  circular- 
ly polarized  wave,  respectively;  the  sign  of  nonlinearity 
is  accounted  by  the  factor  s = sgn(l-p)=  ±1  in  the 
nonlinear  term;  and  k = -rA  12,  rA  = vA  / (%  • 


Eq.  (1)  with  A,  or  A2  is  not  completely  integrable. 
Therefore,  excluding  the  stability  and  asymptotic  analy- 
sis we  used  numerical  integration  for  study  of  evolution 
of  solitons  and  their  collision  dynamics  using  the  special 
simulation  codes. 


2.  Stability  of  2D  and  3D  solutions 


To  study  stability  of  the  GKP  equation  solutions,  we 
performed  coordinate  transformation  and  rewrite  Eqs. 
(1,2)  into  the  Hamiltonian  form 


d,u  = dx(S  H / 8m)  , 


where 


H ^l\-UdxUf+Hd2uJ  +I(v1a;tv)2-«3 


(2) 


Mr, 


d2v  = u , e = piyrl/2,  X = sgn y . The  stationary  solu- 
tions of  Eq.  (2)  are  defined  from  the  variation  problem, 
8 (H  + vPx)  = 0,  where  Px  -jju2dr  is  the  momentum 

projection  onto  the  x axis,  v is  the  Lagrange's  factor, 
illustrating  the  fact  that  all  finite  solutions  of  Eq.  (2)  are 
the  stationary  points  of  the  Hamiltonian  for  fixed  Px  . 
Thus,  conforming  with  Lyapunov's  theorem,  it  is  needed 
to  prove  the  Hamiltonian’s  boundedness  (from  below) 
for  fixed  Px. 

Let's  consider  in  real  vector  space  R the  scale  transfor- 
mations K(x,rx)-»£"l/2r|(1_rf)/2H(;t/£,rx/r|)  (whe- 
re d is  the  problem  dimension,  and  qe  R)  conserving 
the  momentum  projection  Px.  The  Hamiltonian  as  a 
function  of  parameters  t,,  T)  takes  a form 


H (C,iD  = aC~2  +Z>CV2  -<C‘/2t f~d)'Z  +<C4  (3) 

where 

a = -(e/2)j(dxu)2dr,  b=(l/2)j(V±dxv)2dr,  c=ju3dr, 
e = (k/2)j(dxu)2dr . In  2D  case  [d=2  in  expression 
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(3)]  one  can  obtain  that  for  X = 1,  e < 0 the  Hamiltonian 
at  fixed  Px  is  bounded  from  below,  and,  hence,  the  2D 
solitons  are  absolutely  stable.  In  cases  X.  = 1,£>0  and 
X = -1,  e < 0 H has  local  minima,  and  Eq.  (2)  may  have 
the  locally  stable  solutions  for  some  parameters.  All 
other  cases  correspond  to  unstable  2D  solutions. 

In  3D  case  we  obtain  that  the  absolutely  stable  3D  solu- 
tions take  place  for  X - 1,  £ > 0 , and  the  locally  stable 
solutions  can  be  observed  for  X = 1,  £ < 0 if  the  condi- 
tion ab2e/c*  <9/512  is  satisfied.  The  analysis  to  the 
problem  of  the  FMS  waves  beam's  propagation  in  mag- 
netized plasma  enables  us  to  prove  [1],  for  example,  that 
the  3D  beam  propagating  at  0 angle  to  the  magnetic 
field  is  not  focusing  and  therefore  becomes  stationary 

and  stable  within  the  cone  0<arctan  (M  / m)U2  when 

(ml M -cot2  0) 2 [cot 4 0(1  + cot 2 0)]"1  > 4/3 . We  also 
note  that  obtained  results  give  us  the  possibility  to  inter- 
pret correctly  some  numerical  and  theoretical  results  on 
the  dynamics  of  the  internal  gravity  wave  solitons  in- 
duced by  the  pulse-type  sources  in  the  F-region  of  the 
ionosphere  [3]. 

To  study  stability  of  the  3-DNLS  equation  solutions  we 
used  the  formal  change  w -+  h and  investigated  the 
boundedness  of  the  Hamiltonian  [2] 

H = L [ 2 H4  + + jk(V.l3a.vi'J2 ] (4) 

d2w-h,  <p  = arg(/i), 

for  transformation  /i (xr,rx)  — > ^ "i/2ti"i/i(a7^,  r± /r|) 
(£,t|eC)  conserving  Px,  in  complex  vector  space  C. 
The  Hamiltonian  as  a function  of  T)  is  given  by 

H(C,tl)  = aCV2+fcC''  + cCV2  (5) 

where  a =(1  / 2)  Jl /i  I4  dr,  b=Xs  jhh*dx<pdr , 
c=  (a/  2)J (VJL3A.u,)2</r . Solving  the  extremum  problem 

for  functional  (5)  we  obtain  that  Hamiltonian  (4)  is 
bounded  from  below,  i.e., 

H >-3bd/(l  + 2d2),  b< 0 (6) 

if  ac~l  <d  = {l-j2^  -^13  + -Jl&5  , and  in  this  case  3D 
solutions  of  3-DNLS  equation  are  stable.  The  solutions 
are  unstable  in  the  opposite  case,  ac~l  >d,  b<  0.  Con- 
dition b<  0 corresponds  to  the  right  circularly  polar- 
ized wave  with  p = 4mT  I B2  >\ , i.e.  when 
A = l,  s = — 1 , and  to  the  left  circularly  polarized  wave 
when  X = -1,  s = 1 . It  is  necessary  to  note  that  the  sign 
change  A.  = 1 — > — 1,  s = -l— >1  is  equivalent  to  the 
change  /-+-/,  k — » -K  and  for  negative  k the  Hamil- 
tonian becomes  negative  in  the  area  "occupied"  by  the 
3D  wave  weakly  limited  in  the  k L -direction;  in  this 
case  condition  (6)  is  not  satisfied.  Change  of  sign  of  b to 


positive  (when  A = l,  s = l or  X.  = -l,  s = -l)  is 
equivalent  to  the  analytical  extension  of  solution  from 
real  y,  z to  pure  imaginary  values:  y — > -iy,  z — » -iz 

and,  therefore,  equivalent  to  the  change  of  sign  of  k in 
the  basic  equations.  In  this  case  instead  of  inequality  (6) 
the  opposite  inequality  will  take  place.  From  the  physi- 
cal point  of  view  this  means  that  if  the  opposite  ine- 
quality is  satisfied,  the  right  polarized  wave  with  the 
positive  nonlinearity  and  the  left  polarized  wave  with 
the  negative  nonlinearity  are  stable. 

3.  Numerical  results 

In  our  numerical  experiments  we  have  investigated  two 
types  of  the  low-frequency  oscillations  in  a plasma  with 
P«1  and  P>1  which  correspond  to  two  types  of 
waves  and  can  lead  to  the  formation  of  the  multidimen- 
sional solitary  wave  structures.  As  a result,  we  have 
obtained  that  for  FMS  waves  the  2D  and  3D  soliton 
formation  can  be  observed.  In  particular,  we  observed 
the  formation  of  a stable  soliton  with  oscillating  as- 
ymptotics, that  corresponds  to  above  mentioned  analyti- 
cal results.  It  is  interesting  to  note  that  the  2D  soliton 
interaction  dynamics  is  not  trivial  for  GKP  equation 
unlike  usual  KP  equation  [4],  So,  for  example,  for 
X = 1 , E > 0 the  formation  of  a stable  two-soliton 
structure  (so-called  "bisoliton")  can  be  observed  as  a 
final  result  of  interaction  of  two  initial  pulses.  In  the  3D 
case  for  the  FMS  wave  beam  having  the  small  angular 
distribution,  the  stationary  propagation  may  be  observed 
as  a result  of  the  nonlinear  beam  stabilization. 

In  the  case  of  Alfven  waves  propagating  along  the  mag- 
netic field  lines,  we  have  obtained  that  3D  stable  solu- 
tions may  be  observed,  with  3D  spreading  and  collaps- 
ing ones.  These  results  can  be  also  interpreted  in  terms 
of  the  self-focusing  phenomenon  for  the  Alfven  waves' 
beam  as  the  stationary  beam  formation,  scattering,  and 
self-focusing.  Let's  note  that  we  observed  the  dynamics 
of  the  Alfven  waves'  beam  propagating  in  a plasma  with 
p > 1 at  angles  near  0°  with  respect  to  the  magnetic 
field,  and  the  dynamics  of  the  FMS  wave  beam  propa- 
gating in  plasma  with  p«l  at  angles  near  nil  with 
respect  to  the  magnetic  field.  Let's  note  that  for  all  cases 
the  analysis  of  the  Hamiltonian  H deformations  on  the 
numerical  solutions  confirmed  the  stability  of  solutions 
considered  above. 
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